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Making use of the gluon gauge- invariant two-point correla- 
tion function, recently determined by numerical simulation on 
the lattice in the quenched approximation and the stochastic 
vacuum model, we calculate the elementary (parton-parton) 
amplitudes in both impact-parameter and momentum trans- 
fer spaces. The results are compared with those obtained 
from the Kramer and Dosch ansatz for the correlators. Our 
main conclusion is that the divergences in the correlations 
functions suggested by the lattice calculations do not affect 
substantially the elementary amplitudes. Phenomenological 
and semiempirical information presently available on elemen- 
tary amplitudes is also referred to and is critically discussed 
in connection with some theoretical issues. 

13.85.Dz, 12.38.Gc, 12.38.Lg 



I. INTRODUCTION 

The theoretical investigation of quark-quark scatter- 
ing is one of the topical problems in high-energy hadron 
physics. In principle, the possibility to predict this am- 
plitude directly from quantum chromodynamics (QCD) 
could lead to important insight concerning high-energy 
soft processes and meanly elastic hadron scattering. Al- 
though some phenomenological and semiempirical infor- 
mation already exists j^j^ , the determination of this am- 
plitude from a pure QCD approach (model independent) 
is still an open problem. However, recently, remarkable 
progress has been achieved, starting with the nonper- 
turbative approach by Landshoff and Nachtmann in 
which quarks are coupled to Abelian gluons (gluon con- 
densate). In a non- Abelian model, the amplitudes of 
quarks at high energies were computed by Nachtmann 
through the use of the eikonal method Analogous 
version was then developed by Kramer and Dosch jsl in 
the context of the stochastic vacuum model (SVM) 
This model describes the low frequencies contributions in 
the functional integral of QCD in terms of a stochastic 
process by means of a cluster expansion and gives a 
good description of the heavy quark potential The 
most general form of the lowest cluster is the gauge in- 
variant two-point field strength correlator 
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where z = x — y is the two-point distance, a is a char- 
acteristic correlation length, k a constant, g'^{FF) the 
gluon condensate, and Nc the number of colors, C,D = 
1, N^ — l ■ The two scalar functions D and Di describe 
the correlations and are normalized as D{0) = Di(0) = 1. 

These correlation functions play a central role in the 
application of the SVM to high-energy scattering. Once 
one has information about D and Di, the SVM leads to 
the determination of the elementary quark-quark scat- 
tering amplitude and this is the point we are interested 
in. 

Numerical determinations of these correlation func- 
tions, in limited interval of physical distances, exist from 
lattice QCD in the quenched approximation. The first 
determination, by Di Giacomo and Panagopoulos, cov- 
ered the interval 0.4 to 1.0 fm Recently, through 
improved technique (larger lattice) the range was ex- 
tended down to 0.1 fm by Di Giacomo, Meggiolaro, and 
Panagopoulos [ p^ and this brought new and important 
insights on the subject, as will be discussed. 

On the other hand, Kramer and Dosch introduced a 
suitable ansatz for the correlation function D and Di 
which is in agreement with the early lattice results in 
the range 0.4-1.0 fm |10|. Using the SVM they acchieved 
good descriptions of the experimental data on total cross 
sections and slopes of the elastic amplitudes in hadronic 
processes Also, with the same ansatz, Grandel and 

Weise calculated differential cross sections through the 
eikonal approximation (multiple diffraction Glauber the- 
ory). Introducting a monopole parametrization for the 
form factors, a satisfactory description of experimental 
data on pp and pp elastic scattering at small momentum 
transfer was obtained |Q . 

In this paper we make use of the recent lattice results 
for the correlation functions and, through the SVM, we 
calculate the gauge invariant elementary amplitudes in 
both impact-parameter and momentum transfer spaces. 
The results are then compared with those obtained 
through the Kramer-Dosch ansatz and also with some 
phenomenological and semiempirical information avail- 
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able. Since in the SVM the amphtudes are characterized 
by scattering amphtudes for Wilson loops in Minkowski 
space, they will generically be referred to as elementary 
parton-parton amplitudes. 

The material is organized as follows. In Sec. || we 
briefly recall the essential formulas of the SVM connect- 
ing the field strength correlator (P with the parton- 
parton amplitudes and in Sec. |l| we review the corre- 
lation functions from lattice calculations and the ansatz 
by Kramer and Dosch. With this information in Sec. IV 
we present our calculations and results in some detail. 
Discussions and critical remarks concerning nonperturba- 
tive QCD, phenomenological, and semiempirical results 
on parton-parton amplitudes are presented in Sec. ^ and 
conclusions in Sec. VI. 



II. 



PARTON-PARTON AMPLITUDES IN THE 
STOCHASTIC VACUUM MODEL 



In the nonperturbative QCD framework referred to in 
the last section, the study of the elementary scattering 
is based on the amplitudes of quarks moving on light- 
like paths in an external field. In the Nachtmann ap- 
proach the quarks involved in a scattering pick up an 
eikonal phase in traveling through the nonperturbative 
QCD vacuum. In order to have gauge invariant Dirac 
wave function solutions a Wilson loop is proposed to rep- 
resent each quark. In this context the no-color exchange 
arton-parton (loop-loop) amplitude can be written as 



7 = (Tr[7'e"*^^-opi'"""^^'^^"^'"^ - 1] 
X Tr[7'e~^^/-p='"''"^-^^^'"^-l]), 
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where () means the functional integration over the gluon 
fields, the integrations are over the respective loop ar- 
eas, and w is a common reference point from which the 
integrations are performed. 

This expression is simplified in the Kramer and Dosch 
description by taking the Wilson loops on the light-cone. 
In the SVM the leading order contribution to the ampli- 
tude is given by ||] 



7(5) = 7,62(6), 



(3) 



where b is the impact parameter, 77 is a constant depend- 
ing on normalizations [see Eqs. (|22) and (p3)], and 



J J daf^^dap^TT{Ff_,^{x;w)Fp^{y;w)). (4) 

Here {g'^F^^{x; w)Fp„{y; w)) is the Minkowski version of 
the gluon correlator. We will return to this point later. 

After a two-dimensional integration, eih) may be ex- 
pressed in terms of the correlation functions in Eq. (|^) 

by [| 



e(&) = ei(6)-f eii(6), (5) 
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where is the four-momentum transfer squared and D = 
D or Di, 



(8) 



with Tn denoting a n-dimensional Fourier transform. 

In the impact parameter space 7 represents the elemen- 
tary profile function, from which the elementary parton- 
parton scattering amplitude is calculated through a two- 
dimensional Fourier transform: 



f{q') = / bdbMqbhib), 
Jo 



(9) 



where Jq is the Bessel function. 

With the above formalism, once one has inputs for the 
correlation functions D{z) and Di{z) the elementary am- 
plitudes in the impact parameter and transfer momentum 
spaces, Eqs. (|^) and (||), respectively, may, in principle, 
be evaluated through Eqs. (p[)-(p[). 



III. CORRELATION FUNCTIONS 

We now recall the theoretical information available 
about the correlation functions D and _Di, namely, that 
obtained from lattice QCD and the ansatz introduced by 
Kramer and Dosch. In S 

in order to calculate the elementary amplitudes. 



IV] we use these functions 
The 

translation from Euclidean to Minkowski space will be 
discussed in the next section. 



A. Lattice results 

The first determination of the correlation functions 
from lattice QCD in the quenched approximation was 
obtained through the cooling technique with a 

lattice of size 16^ and in the interval of physical distance 
0.4 to 1.0 fm ||ll[. In this range the theoretical data 
showed agreement with an exponential decrease of both 
correlators with the distance. In particular, using the 
notation of Ref. pO| , the function D is parametrized by 
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where r is the physical distance (Euclidean space), and 
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with a typical statistical error of a few percent. In Ref. 
PI the value of A was determined from string tension 
studies M| leading to 



C = 137 fm" 



A = 0.22 fm, 



(12) 



which corresponds to A = 4.9 MeV. 

Recently new results were obtained on a 32"* lattice 
allowing the determination of the correlators at distances 
down to 0.1 fm jl^. A novel result was the appearance 
of a deviation from the exponential behavior at shorter 
distances (< 0.4 fm), indicating a 1/x'* divergence at the 
origin. Putting together the data of Refs. [Ill and [ p2[ , 
good agreement is obtained within statistical errors, with 
a parametrization of the form [[l2| 
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As before, from string tension information on A, Xa = 
0.22 fm and As = 0.43 fm. The typical error is again of 
the order of a few percent. In Sec. IV we test and discuss 
the influence of the divergent terms. 



B. Kramer-Dosch ansatz 

In the approach of Ref. ||l^, the effect of _Di(z^/a^) is 
neglected since, from QCD lattice calculation, the value 
of K in Eq. (|^) is « 1 ||ll|. So, we will consider only D 
for this case. 

Assuming functions which have a well defined Fourier 
transform and that can be analitically continuated to 
Euclidean world in the momentum space, Kramer and 
Dosch introduced a family of functions depending on an 
integer parameter n . These functions may be analiti- 
cally integrated and the selection of the n value was made 
by comparison with the early lattice results, Eq. (10). 
In the Euclidean space the fit in the interval 0.5-0.8 fm 
(early lattice results) led to n = 4 and for A = 4.4 MeV 
it was obtained a = 0.35 fm and K{g^FF) = 1.774 GeV* 
The final result for the correlation function with 
n = 4 is given by 
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where Ki , Kq are Bessel functions and x is a dimension- 
less variable 

37r z 



C. Remarks on the correlation functions 

In what follows, it is important to stress that the 
Kramer-Dosch ansatz for the correlation function D is 
an analytical form, firstly selected by analytical continu- 
ation and asymptotic limit conditions and then through 
fit to early lattice information, Eq. (|lo|). On the other 
hand, "real" lattice results correspond to sets of discret 
theoretical points with errors, in a finite interval of phys- 
ical distances (0.1 - 1.0 fm). The parametrizations intro- 
duced in both early ||ll| and recent jlj] results, extrapo- 
late this interval down and above. Although, presently, 
there is no theoretical informations on these extrapo- 
lated regions we will assume the parametrizations ( p^ ) 
and (^^ as representing lattice results and will refer as 
lattic e para met rizati ons. We will return to this point in 
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The results reviewed in this section for the correlation 
function D are displayed in Fig. |^ for comparison and this 
deserves some discussion. We plotted the curves corre- 
sponding to the early lattice result, Eqs. ( [lO| ) and (|T^), 
the new lattice result, Eq. (13), and the Kramer-Dosch 
ansatz, Eq. (p^. As mentioned above, the final form of 
this ansatz was obtained through a fit to the early lattice 
parametrization in the interval 0.4 — 0.8 fm. From Fig. 
|l|, there is a discrepancy between these two curves in the 
above interval, which does not appear in Fig. 10 of Ref. 
Ip^ . The reason is the different values of the parameter 
A, namely, 4.4 MeV in Ref. and 4.9 MeV in Ref. 
Ill and Fig. This difference is significant since, for 
example, the parameter C in Eq. (|ll|) is proportional to 
A''. On the other hand, the differences between early and 
new lattice results in the above interval are compatible 
with the errors associated with each curve. 



IV. PARTON-PARTON AMPLITUDES 

In the last section we reviewed the theoretical informa- 
tion presently available concerning the correlation func- 
tions, namely, the recent lattice parametrizations, Eqs. 
( p^ ) and (p^, and the result from the ansatz by Kramer 
and Dosch, Eq. (|l5|). 

In the context of the SVM (Sec. Q), the corre- 
sponding predictions for the parton-parton amplitudes 
are obtained, in principle, through four steps, envolv- 
ing the following calculations: (1) D{k) and Di(k) by 
the four-dimensional Fourier transform, Eq. (^); (2) 
two-dimensional inverse Fourier transforms of D{k) and 
dDi/dk^, which enter in Eqs. (^) and (0); (3) eikonal 
phase through Eqs. (^ and (^; (4) elementary ampli- 
tudes in the impact parameter space (profile function), 
Eq. (^, and in the transfer momentum space (scatter- 
ing amplitude), Eq. (^). These are the central points of 
this work and in this section we present and discuss the 
calculations in some detail. 
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FIG. 1. Correlation function D given by (a) recent lat- 
tice parametrization, Eq. (jl3|) (dotted); (b) early lattice 
parametrization, Eq. (|l^) (dashed); (c) Kramer and Dosch 
ansatz, Eq. ( p^ (dot-dashed). 



A. Results from lattice parametrizations 

Our basic assumptions in using lattice results are the 
following. 

(a) We consider the parametrizations ( |l3| ) and (|lj) as 
representing the correlation functions and so, the corre- 
lation lenght has the value a = 0.22 fm. 

(b) The effect of the extrapolation down to 0.1 fm will 
be investigated by taking account or not of the divergent 
term 1/|2:|^, that is, taking or not B = and Bi — in 
Eqs. ( p^ ) and (14), respectivelly. 

(c) We consider the lattice parametrizations, Eqs. dl^ ) 
and (^), as representing the correlator functions in the 
Euclidean world. With this assumption the above corre- 
lation functions enter into Eq. (|l]), with the adequated 
tensor structures, and then directly in Eq. (|8|). 

In using the above assumptions we seek to see the con- 
sequences in terms of the elementary amplitudes. 

The first step referred to before concerns the four- 
dimensional transform 



V{k) = F4[V{z)] = J d^zViz) cxp(ik • z) 



(16) 



for T) = D, Di, here in Euclidean space according to 
assumption (c) above. Although the lattice data arc lim- 
ited to the interval 0.1-1.0 fm, the parametrizations ( p^ ) 
and (|lj) extend to all the space and include the divergent 
term l^;!"**. However, we found that when the lower inte- 
gral limit Zm becomes smaller than w 10"'^ fm, numeri- 



cal evaluation of this transform , for both correlators, 
may be put in the form 



Vik)=d{k)+C{z,n), 



(17) 



where d{k) — di{k), d{k) are smooth (finite) decreasing 
functions of fc, and C{zm) is a constant which increases 
when the lower integral limit becomes smaller. This effect 
is shown in Fig. ||. 

5.5 
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FIG. 2. Four-dimensional Fourier transform of D{z), Eq. 
(p^, and the effect of the lower integral limit 2,„ (between 
10"^ and 10"'^) as represented by Eq. (0): (a) d{k) = d(fc); 
(b) d{k) = di{k). 
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Now, in the second step, we should calculate 
the two-dimensional Fourier transforms of D{k) and 
dDi{k'^) / dk"^ . Because of the derivative in the last case 
the constant C{zm) may be neglected. In the former case, 
since the Fourier transform of C{zm) leads to a (5 func- 
tion, its effect in the calculation of the eikonal phase e/ 
in Eq. (||) has also no influence. 

With this, earring out the numerical integration down 
to 10~^ fm we have a stable behavior for d{k) in both 
cases (see Fig. ||). The numerical results for d{k) and 
di{k) are shown in Fig. |[ We then proceed to fit these 
points through the CERN-MINUIT routine by func- 
tions of the type 

2 

d{k) = aj exp{—bjk) + as exp(— 63^:^), 
2 

di{k) = ^aij exp(-6ij/c) -|- 013 exp(-6i3A:^), (18) 

also shown in Fig. ||. The values of the free parameters 
are displayed in Table |. 
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FIG. 3. Finite part of the correlation functions in momen- 
tum space, Eg. (Ilj), calculated numerically (points) and fits 
through Eq. ( |l^ ) (solid curves). 



TABLE I. Values of the free parameters for d{k) and di(fc) 
in Eq. (p^. Parameters fei, 62 and bn, bi2 are in fm while 
63, bis in fm^. 



j 






ai,/(4^2) 


bij 


1 


0.50003 


0.094168 


0.46314 


0.12874 


2 


0.76546 


0.010859 


0.28780 


0.48787 X 10"^ 


3 


0.73263 


0.104490 


0.18440 


0.36088 X 10"^ 



With Eq. ( p^ and the derivative of di{k), step 2 is 
performed analytically leading to 

2 

^2'[d{-q^m) = 5]a,[/3,62 ^ ip3/2 _^ g^p(_^^^2)^ 

(19) 



(6)=^ai,[/?y62 + l]-i/2 + ^^3 

xexp(-/3i362). (20) 

From Eq. <^^), in step 3, the contribution ei{b) to the 
eikonal phase is analytically calculated: 



ei(6) = Ciexp(-ei6') + C2[Erfc(e2fe)]& + 5]a,([(e,6)2 

+ l]i/2 + £,6), (21) 

where Erfc is the complementary Error function. The 
contribution en (6), Eq. (^, is given directly by Eq. ([2^). 
Figure ^ shows both ei and en as function of the impact 
parameter. We will discuss these results in Sec. 

At last, step 4 is performed with e(6) = [ei(&) + eii(&)] 
into Eq. (||), leading to the elementary profile 7(&), Eq. 
(0), with normalization constant 



'n 



9 X 82' 



(22) 



Then numerical integration furnishes the corresponding 
elementary parton-parton amplitude /((?), Eq. (9|). The 
results for 7(6) and /(q) are shown in Figs. |5 and |^, 
respectively. 

As expressed in the beginning of this section, in or- 
der to investigate the effect of the divergent term in the 
recent lattice results, we also calculated 7(6) and j{q) 
neglecting these terms in both correlators, i.e., by taking 
account only of the first term in Eq. ([l3| ) and (|lj). All 
these results (with and without the divergent terms) are 
displayed in Figs. || and|| for 7(6) and /(q), respectively. 
We will return to this point in Sec. ^ after the discussion 
of 7(6) and /(g) predicted from the Kramer and Dosch 
ansatz. 



^For short, in what follows, we will refer only to the analyt- 
ical structure of the parametrized and calculated functions. 
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FIG. 4. Contributions for the eikonal phase, Eqs. (^), 
and (0), from lattice (ei and en from D and Di, respectively) 
and Kramer-Dosch ansatz (ekd)- 



B. Results from Kramer-Dosch correlators 

Concerning step two, a suitable parametrization for 
the two-dimensional inverse Fourier transform of D{k) 
was introduced in Ref. ]lO| , so that step three may be 
performed analytically leading to 



ckd = IC exp(-x) ^ QnX''' 

^2 pz?\„4ol4 /fo4:^3\ 



(23) 



where /C = K{n^FF)a'^2^'^/{3^Tr^). This result is dis- 
played in Fig. ^ together with those obtained through 
the recent lattice results. 

At last, as before, the elementary profile is calculated 
through Eq. (0) now with normalization constant 



1 



9(8 X 12) 



(24) 



By numerical integration we obtain the scattering 
amplitude. The results are shown in Figs. ^ and|6| 
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FIG. 5. Elementary (parton-parton) profile functions, Eq. 
(^): (a) from lattice parametrizations with the divergent term 
(dotted), without the divergent term (dashed); (b) from the 
Kramer-Dosch ansatz. 



V. DISCUSSION 



A. Partial conclusions 



In this section we first sketch the conclusions coming 
directly from our calculations and then present discus- 
sions related to some phenomenological/semiempirical 
informations available on elementary amplitudes and also 
some theoretical issues. 



The main results of our calculations are displayed in 
Figs. 1^, 1^ and || (corresponding to eikonals, profiles, and 
amplitudes, respectively) and lead to the following con- 
clusions. 
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tions, there is no agreement at all between profiles or 
amplitudes from lattice parametrizations and from the 
Kramer-Dosh ansatz. 
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FIG. 6. Normalized elementary (parton-parton) scattering 
amplitudes, Eq. (^), from lattice parametrizations with the 
divergent term (dotted), without the divergent term (dashed), 
and from the Kramer-Dosch ansatz (dot-dashed). 

From Fig. ^, the eikonal phase e — ei + en from 
the recent lattice parametrizations (including the diver- 
gent term), has a maximum at 6 = and decreases 
smoothly as the impact parameter increases. Differently, 
the eikonal phase from the Kramer-Dosch ansatz presents 
the maximum at 6 « 0.5 fm and reaches zero at 6 = 0. 
These distinct behaviors come from the differences be- 
tween the correlators, shown in Fig. |l]. 

Due to the structure of Eq. (||), the corresponding 
profiles present similar behaviors, as shown in Fig. ^. In 
this case we investigate the effect of the divergent term, 
referred to before. We see that, with or without the diver- 
gent term l/l^l'* in the correlators D and Di, the profiles 
present maxima at 6 = and a smooth decrease with the 
impact parameter. With and without the divergence, the 
profiles reach ~ 10% of there maxima at 6 = 0.1 and 0.2 
fm, respectively. We conclude that the term 1/|2;|'' does 
not significantly alter the outputs leading only to a more 
concentrated and smaller profile near 6 = 0. 

As a consequence of the profiles, the corresponding 
scattering amplitudes also present analogous similarities 
and differences, as can be seen in Fig. ^. Normalization of 
the amplitudes to 1 a±q^ = 0, shows that both lattice am- 
plitudes (with and without the divergent term), present 
a similar smooth decrease as the momentum transfer in- 
creases, going to zero through positive values. The effect 
of the divergent term is a faster decrease only. Differ- 
ently, the amplitude from the Kramer-Dosch ansatz de- 
creases still faster, presenting a change of sign (zero) at 
g2 ~ 0.5 GeV^ and goes asymptotically to zero through 
negative values. We see that according to our calcula- 



B. Phenomenological and semiempirical 
informations 

Presently, some limited informations concerning el- 
ementary profiles and/or amplitudes are available 
from phenomelogical models and semiempirical analysis 
(model independent) as explained in what follows. The 
usual framework of analysis is the eikonal approximation 
in which the elastic hadronic amplitude F is related to 
the eikonal function x by |l^ 

Fiq, s)^i j bdbJoiqb) {1 - exp[ix(fo, s)]} , (25) 

where ^/s is the center-of-mass energy. The hadronic am- 
plitude is connected with physical observables as differ- 
ential, total cross sections, etc., and the eikonal is defined 
in the context of each phenomenological model. In the 
simplest approach, represented by the Glauber or Chou- 
Yang formalism, the eikonal is expressed by [OlPO|j21 



X{b,s)^C / qdqJoiqb)GAGBf, 



(26) 



where Ga, Gb are hadronic form factors, / the elemen- 
tary (parton-parton) amplitude, and C depends only on 
the energy. 

Phenomenology. In the absence of theoretical pre- 
dictions for both form factors and elementary ampli- 
tudes, models are yet characterized by different choices 
of parametrizations for these functions. As reviewed in 
Ref. ||l|, some models assume parametrizations for the 
electromagnetic form factors and others introduce energy 
-dependent form factors. With quite different choices for 
the elementary amplitudes all these models present sat- 
isfactory descriptions of the experimental data. For this 
reason it is difficult to identify a "phenomenological am- 
plitude" since it depends on the parametrizations used 
for the form factors and some other specific aspects of 
each model. So this puts serious limitations in attempts 
to compare theoretical - phenomenological results and we 
will return to this point in the next section. 

However, even taking account of these limitations, it 
may be useful to see what kind of similarities could be 
estabilished. In this sense, based on the review of Ref. 
||l| we could say the following. 

(a) The behavior of both 7(6) and f(q) predicted from 
lattice parametrizations (Figs. |^ and P) are in qualita- 
tive agreement with the form introduced by Glauber and 
Velasco I2I 



/Gv(g,s) 



1 



r 



(27) 
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where c? = 7.14 GeV^ for "pp scattering at -/s 
GeV (with a phase factor in /gv) and c? = 0.625 GeV 
for at 546 GeV. The elementary amphtude decreases 
smoothly through positive values and goes asymptoti- 
cally to zero, as shown in Fig. for comparison with the 
results from lattice parametrizations. 
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FIG. 7. Elementary amplitudes from lattice parametriza- 
tions, with (dotted) and without (dashed) the divergent terms 
(same as Fig. ^), cornpared with the Glauber and Velasco 
parametrization, Eq. (^7|), for pp at 23 GeV (solid) and pp at 
546 GeV (dot-dashed). 

(b) The profiles and amplitudes from the Kramer and 
Dosch ansatz (Figs. || and g) show qualitative agreement 
with the parametrizations used by Menon and Pimentel 
i3: 



niBSW — 



1 - q^ja^ 



(28) 



where o? = 8.2 GeV^, independent of the energy or reac- 
tion, and the subscript means modified BSW (Bourrely, 
Soffer, and Wu), as explained in Ref. The amplitude 
presents a change of sign and then goes asymptotically 
to zero through negative values. 

Semiempirical analyses. This last behavior is also pre- 
dicted in the semiempirical analysis by Buenerd, Furget, 
and Valin [|j . These authors used the fit to pp differential 
cross section data introduced by Amaldi and Schubert 
and calculated the eikonal in the transfer momen- 
tum space. In the context of the Glauber-Chou-Yang 
model, Eq. (26), this quantity reads 



xiq,s)^CGlf. (29) 
The elementary amplitude /(g) was then extracted at 



■\/s = 23.5 GeV by dividing x(g, s) by parametrizations 
for the electromagnetic form factor (Feist parametriza- 
tion and a dipole form used in the Bourrely-Soffer-Wu 
model 1^) and normalizing the result at g^ = to one. 
The amplitude so extracted presents a change of sign at 
Qq = 8.6 GeV^, going to zero through negative values. In 
the CERN Intersecting Storage Ring (ISR) energy region 
(23.5 GeV < < 62.5 GeV) the position of the zero de- 
creases as the energy increases, reaching — 5.0 GeV^ 
at yi= 62.5 GeV |§. 

However, this kind of procedure has two critical as- 
pects. 

(1) Since experimental data are available only in a lim- 
ited interval of momentum transfer (q^ < 6.0 GcV^ in 
the above case) all kinds of extrapolations in the fits, al- 
lowed statistically, should be taken into account. This 
can be made by error propagation from the fits parame- 
ters, which, however, was not done in the above analysis. 

(2) The result for the amplitude yet depends on 
parametrizations for the form factor. The choice to use 
an electromagnetic form factor is an approximation since 
we are treating a hadronic and not an electromagnetic 
interaction. However, there is no theoretical or experi- 
mental information on hadronic matter form factors. 

Recently, through fit procedures to pp differential cross 
section data, Carvalho and Menon obtained the eikonal 
in the transfer momentum space, taking account of error 
propagation and also the effect of large momentum data 
p5[ . The result shows statistical evidence for a change of 
sign in the eikonal in the interval 5 GeV^ ^ 9^ ^ 9 GeV^ 
at the ISR energy region |2^ and asymptotical limit to 
zero through negative values. However, the movement 
of the zero with the energy cannot be inferred on sta- 
tistical grounds. The analysis was performed through 
both a fit method and a numerical method. The for- 
mer approach gives an average position of the zero at 
ql_= 7.1 ± 0.7 GeV^ and the later at g^ = 6.1 ± 0.7 GeV^ 



5|. If, in the context of Eq. (|29|), this zero is not asso- 
ciated with an hadronic form factor the result favors the 
amplitude from the Kramer-Dosch ansatz. Moreover, we 
found that the above positions of the zero are well repro- 
duced if the gluonic correlation length in the Kramer- 
Dosch amplitude has the value a ~ 0.1 fm. Specificaly, 
a = 0.0957 fm for ql = 7.1 GeV^ and a ^ 0.1040 fm 
for (7q — 6.1 GeV^. The result, in the case of the zero 
at gg = 7.0 GeV^, is shown in Fig. ||, together with 
the parametrization used by Menon and Pimentel for 



= 7.0 GeV^ in Eq. 



C. Theoretical and phenomenological issues 

The physical picture in the Nachtmann approach is 
that of interaction of two long-lifetime partons with 
a common nonperturbative QCD vacuum, the partons 
traveling essentially on straight lightlike world lines. It is 
also assumed that the transversal momenta of the partons 
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(related to the hadron's movement) are orders of magni- 
tude smaller then the forward momenta. Therefore, this 
picture is characterized by very small momentum transfer 
and asymptotically high- energy regime. These are crucial 
aspects of the nonperturbative approach: They specify 
detailed frontiers, beyond which physical interpretations 
and/or comparisons with phenomenological information 
have no justification and this demands further discussion. 
In Figs. 1^ and ||, we display the results for the elementary 
profiles and amplitudes in a wide interval of impact pa- 
rameter and momentum transfer. Based on the range of 
validity of the theoretical framework, we must see part of 
these curves as extrapolations to large momentum trans- 
fer and an essentially small impact parameter. Also, the 
elementary profiles and amplitudes do not depend on the 
energy since they are associated with the asymptotically 
high-energy limit. This leads to the following considera- 
tions. 
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FIG. 8. Elementary amplitudes from Kramer-Dosch ansatz 
with correlation lenght a = 0.0957 fm (dot-dashed) compared 
with the parametrization used by Menon and Pimentel, Eq. 
(^, with = 7.0 GeV^ (dotted). 

For the amplitudes obtained from the lattice 
parametrizations, the monotonical decrease to zero 
through positive values is clearly an extrapolation, at 
large transfer momenta. This behavior, however, could 
be broken, for example, by a specific model for that re- 
gion, introducting scale properties, and even zeros for 
large . 

The comparison displayed in Fig. ^ between the 
lattice inspired amplitude and the Glauber- Velasco 
parametrization has also only a qualitative character, 
since in the last case the amplitude depends on the en- 
ergy (and reaction). We see, however, a tendency to 



agreement as the energy increases, in the small transfer 
momentum region. 

The general agreement between the behaviors of the 
Kramer-Dosch amplitude and /mBSW (Fig- ^) niay have 
some more profound (but limited) meaning, since both 
amplitudes do not depend on the energy. Also, the con- 
nections with exp erimental data in the approaches of 
Refs. pofl and [p3| are obtained with energy dependences 
in the hadron's transverse wave function and hadronic 
form factors, respectively. This is very similar from a for- 
mal point of view: The energy dependence is ultimately 
associated with the hadronic radius. However, the lim- 
its in the transfer momentum region must be taken into 
account. 

Another theoretical limitation that should be stressed 
concerns the correlation functions determined on the lat- 
tice. As referred to in Sec. IV A the theoretical results 
were obtained in the finite interval 0.1-1.0 fm and in this 
paper we used the parametrizations introduced in Ref. 
Q . With this assumption we found that the divergent 
term 1/|2;|'' does not affect the profile or amplitudes in 
a substantial way (Figs. |^, H). Now, what kind of infor- 
mation can we extract from the "unknown" region 0-0.1 
fm? We observed that, neglecting the divergent term in 
the correlation function D{z) by taking i? = in Eq. 
(|l3|), leads to a similar result as that obtained in the 
early calculation, Eq. (|lO|). Specifically, in the former 
case the finite D value at \z\ — is « 4.7 GeV* a nd in 
the latter case, « 5.0 GeV . From Fig. |l] and Sec. |IIIB 
we see that these values are yet near three times the fi- 
nite value of the Kramer-Dosch correlator (w 1.8 GeV^). 
This difference between the lattice extrapolations (with- 
out divergence) and Kramer-Dosch ansatz in the unaces- 
sible region of small distances (0.0-0.1 fm) has strong 
consequences, leading to the differences observed in the 
corresponding profiles (Fig. ||) and amplitudes (Fig. ||). 
That is, the existence of zero in the amplitude seems to 
be associated with the particular form of the correlator 
at < \z\ < 0.1 fm (for example, exponential, Gaussian, 
etc.), or even at 0.1 < \z\ < 0.4 fm, and also to the value 
of the maximum at \z\ — 0. This small distance region, 
however, is not presently acessible to lattice calculations. 

The results for the profiles and amplitudes showed a 
weaker dependence on the divergence of the type \z\~'^ 
than should be expected. Theoretically this kind of diver- 
gence may be associated either to operator product ex- 
pansion (OPE) divergences or some specific high-energy 
gluon exchange process. 

Our comparative analysis also brings some informa- 
tion on the gluonic correlation lenght a that should be 
discussed. This quantity is a very important parameter 
in the SVM applied to high-energy scattering. All the 
formalism is strongly sensitive to its value since, for ex- 
ample, total cross sections for parton-parton scatte ring 
is proportional to a^" [^. As explained in Sec. IIIB, the 



Kramer-Dosch (KD) correlator function was determined 
by fit to early lattice results, from which they obtained 
a = 0.35 fm and this value may reproduce parton-parton 



9 



total cross sections of the order of 4 mb ||]. On the 
other hand, from SU(2) lattice gauge theory, it was esti- 
mated a = 0.1 ~ 0.2 fm |2^], and both early and recent 
parametrizations of the correlation functions from the 
lattice, Eqs. ® , (|l|), and (|l|), lead to a = 0.22 fm. In 
the context of the SVM, a ~ 0.2 fm means a very small 
contribution to parton-parton scattering in case of using 
KD ansatz, as discussed in ||]. 

Concerning these distinct estimates, we showed in 
our analysis that for a = 0.35 fm the elementary am- 
plitude for the Kramer-Dosch correlator has a zero at 
~ 0.5 GeV^ and that for a ~ 0.1 fm the position 
changes to — 7.0 GeV^. Even taking account of 
the limitations related to the semiempirical analysis (Sec. 
VB ), it is evident the existence of zeros at the ISR en- 
ergy region, m the interval 5.0 < < 9.0 GeV^ ||,||. 
We conclude that, although this interval corresponds to 
finite energies, it seems to favor smaller values for the 
correlation lenght than 0.35 fm. 

At last it should also be stressed that the lattice re- 
sults we used were obtained in the quenched approxima- 
tion and therefore additional uncertainties may be con- 
sidered. Recently, dynamical fermions have been taken 
into account in lattice calculations, allowing a full QCD 
treatment p^ . However, since a reanalysis is in course, 
with improved statistics [g9[, we will leave this subject for 
future consideration. It will also be suitable to compare 
these new results with those presented here. 



VI. FINAL CONCLUSIONS 

In the limit of extremely high energies (quarks on the 
light-cone) and small momentum transfer between par- 
tons, the nonperturbative approach connects the correla- 
tion functions with the parton-parton amplitudes. Mak- 
ing use of recent lattice parametrizations and also of the 
Kramer-Dosch ansatz for these correlators, we calculated 
the corresponding elementary amplitudes in both impact 
parameter and transfer momentum spaces. In the case of 
lattice inputs the calculations were performed with and 
without the divergent term recently introduced. 

We found that this divergence in the correlator does 
not substantially affect the normalized amplitudes and 
also that the results from lattice parametrizations and 
from the Kramer-Dosch ansatz do not agree. Although 
the lattice parametrizations are statistically consistent 
with the theoretical points (x^/Ndf 1-7 [Q) there 
is no real theoretical information in the interval 0.0-0.1 
fm. This is a crucial point and our analysis furnished two 
novel results concerning possible connections between the 
short range behavior of the correlators and the final re- 
sults for the elementary amplitudes. On the one hand. 
Fig. ^ shows that the divergent term l/lz]"* causes only a 
faster decrease of the amplitude as the momentum trans- 
fer increases, that is, it does not generate local minima 
or zeros. The amplitude coming from this particular lat- 



tice parametrization, with or without the divergent term, 
is characterized by a monotonical decrease with the mo- 
mentum transfer through positive values. On the other 
hand, in the case where the correlator may be finite at 
\z\ = 0, its maximum value and also its shape at small 
\z\ {\z\ < 0.4 fm) strongly affect the behavior of the am- 
plitude, as sh own in Fig. || and quantitatively discussed 
in Sec. |V C| : Zeros (change of sign) and local minima 
may be generated. That is, despite the present lack of 
theoretical information in the short range 0.0 — 0.1 fm, 
we can infer, in the context of the stochastic vacuum 
model, what kind of behavior should be expected for the 
elementary amplitude. 

Attempts to compare these results with phenomenolog- 
ical and/or semiempirical information on the amplitudes 
are limited. The reason is that this information comes 
from an analysis at finite energies (presently available) 
and the amplitudes depend on particular model assump- 
tions or specific parametrizations for the form factors. 

On the other hand, as referred to in Sec. 0, the distinct 
results for the amplitudes from the lattice and Kramer- 
Dosch ansatz show qualitative agreement with different 
model predictions such as Glauber- Velasco and Menon- 
Pimentel, respectively. Furthermore, semiempirical anal- 
ysis could favor the Kramer-Dosch amplitude, but with 
a gluonic correlator length smaller than 0.35 fm. 

We conclude that, eve n tak ing account of all the limita- 
tions referred to (Sees. VB and |V C| ), these agreements 
suggest possible links between theory and phenomeno- 
logical or semiempirical analyses. Further investigations 
along these lines may be important as a source of feed- 
back for theoretical (nonperturbative) ideas, meanly con- 
cerning suitable considerations on energy dependences 
and small distance phenomena. 
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